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CN . We discuss the color-superconductivity and its effect on the cooling behavior of strange 

■ quark stars. The neutrino emissivity and specific heat of quark matter are calculated within 

the BCS theory. In the superconducting phase, the emissivity decreases and causes suppres- 
sion of the cooling rate. It is shown that the phase transition leads to a sudden discontinuous 
. suppression of the cooling rate in cooperation with the specific heat. 

§1. Introduction 

Recently many authors have investigated the color-superconductivity in quark 
matter. EP Although there is no evidence of its existence, it is pointed out that there 
might be exist in the core of neutron stars and in the central region of high-energy 
heavy-ion collisions. 

What physical quantities must we observe in order to show the existence of 
color-superconductors? In this paper, we discuss physical phenomena, which are 
characteristic of the color-superconductivity. To this end^we will focus our attention 
on the cooling behavior of the hot strange quark stars, a' The recent works show 
that the quark matter is transformed from Fermi gas into superconducting fluid with 
the decrease of the temperature; there occurs a phase transition of second order at a 
critical temperature. & Therefore it is expected that some physical quantities have 
^ ■ singularity when the hot quark matter is cooled down. This singularity may be a 

possible signal for the existence of color-superconductors. 

The coolinff phenomenon of quark stars has been studied for many years by 
several authors. t§ According to their considerations, only the week interaction can 
release the energy from such a high-density quark matter; the simplest possible 
process is neutrino emission from the quark matter and in particular the pair of d 
quark (3 decay reaction is dominant (URCA): 

| d —* u + e~ + v e , , , 

I u + e~ — > d + v e . 

The neutrino energy loss at an interval dt for via this process is denoted by C m dt. 
The constant C m is called neutrino emissivity. On the other hand, noting that it 
is expressed by CydT with the specific heat Cy of the quark matter, we obtain a 
differential equation determining the temperature as a function of time, 

dT C m (T) 
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dt C V {T) 
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This equation describes the cooling behavior of the quark matter. 

Recently much efforts have been made in order to study the influence of the 
color-superconductivity on the cooling behavior. B 1 They are, however, restricted 
to the low temperature regions and dealing with mainly the specific heat. It is 
the purpose of this paper to investigate the effect of the neutrino emissivity on the 
cooling behavior around the critical temperature. 

In the next section, the mean field (BCS) theory for the color-superconductivity 
is developed in such a way that it is suitable for our present case of finite temperature 
mean field theory. Then we calculate the neutrino emissivity and specific heat of 
quark matter in §3. Finally (§4) numerical calculations are carried out for a typical 
quark matter and some comments are given. 



§2. BCS theory at finite temperature 



First let us develop the BCS theory in the three-flavor quark system at finite 
temperature. The Lagrangian is assumed to be given by 

C = W(i-y ■ 8 + //7°)!F + ^(^A^X^A^), (3) 

a 

where [i is a chemical potential and A a denotes the color SU(3) matrix. This effec- 
tive Lagrangian comes from the one-gluon exchange interaction with infinite gluon 
mass due to the many-body effect in the medium. To study pair correlation, it is 
convenient to use the Fierz transformation tJ), 



2 w,^- 



C = #(* 7 • d + /i 7 )^ + ^g^2' (W^ 5 CX a A b W t )(^ t C- 1 j 5 X a A b W), (4) 
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where A b denotes the flavor SU(3) matrix. Here we have left only the most attrac- 
tive terms: spin-singlet, color-antisymmetric and flavor-antisymmetric (a, b = 2,5, 7) 
terms. 

The partition function of our system is given by the following functional integral, 

Z = Trexp(-/3i?) = / VWVW exp(- J d 4 xC), (5) 

with the use of (5 = 1/T and d x = drd 3 x. Here let us introduce auxiliary field f{x) 
in terms of the following identity: 



1 = J V<p*Vipexp(-K 2 J d 4 x\ip p {x)\ 2 ). (6) 



Using this equation, we have chosen the constant k in the above function in such a 
way that the four-Fermi interaction term is eliminated, 

Z = J VWV&Vp*Vipexp[- J d 4 x{^(^7 • d + fj,-y )^ 

+(Wj 5 CX a A b W t tp p + $ t C~ 1 'f i \ a A h ¥(p* p ) + k 2 |^| 2 }], (7) 
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where the color-flavor matrix T p is defined by A a ® A b . From the above equation it 
is clear that the auxiliary field is equivalent to a wave function of Cooper pairs. 

Here we take a mean field approximation (BCS theory), assuming that the field 
(p(x) is a constant, 

tpp(x) = 5 a ^h,2m- (8) 

This expression means that our Cooper pairs are u-d spin-singlet with the momenta 
±fc. The reason that we have taken u-d pairing is that the mass of the s quark is 
heavier than those of u,d quarks. 

To carry out the integrations, we expand the Fermi field & into Fourier series: 

^( x ) = /Z Z~2 b p,MP, s)exp{i(p ■ x - uj h t)}, (9) 

n p s 

where the Matsubara frequency oj n represents = (2n + 1)/3 _1 and s (±1) is an eigen- 
value of the quark helicity. Substituting this equation into the partition function, it 
can be represented as the following matrix form: 



(A U \ iuJ n~(p A) W b p l \ 2, ,2 

b -^){ ai ^ + cj^; ) _ 1 ) +K ^ 

(10) 

where the usual gap parameter has been defined by Aq = —2ipoT 2 (T 2 = X 2 ® A 2 ). 
This Cooper pair is antisymmetric in spin, color and flavor spaces independently. 
Carrying out the functional integrals of the quark fields, we obtain 



^TrOog/?^-^))-/^ 



Z = Nexp 

where E p is the energy matrix of the quasi-particle and given by 



4 



(11) 
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E p " = Q 2 V + A 2 {\ 2 ® A 2 ) 2 . (12) 

The characteristic of the quasi-particle is the following; all the s quarks are gapless 
and have not pair correlation, on the other hand the u, d quarks with the color indices 
1 and 2 have a finite gap and the remaining particle is gapless.^ 

From the partition function (7), one can easily calculate the thermodynamical 
potential, 

n = -/r 1 log z = -4^(^ p + 2/r 1 i og (i + exp (-/3E p )) 

p 

-SzZttp + 2/T 1 log(l + exp(-/^)) + -k 2 A 2 . (13) 
p 

The gap parameter can be determined through the stationary condition of the ther- 
modynamical potential, 
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dA~o~ ~ 



4£-^(l-2n p ) + -A) = 0, (14) 
P 1 



with n p = (1 + e^ p ) -1 , the distribution function of the quasi-particle. 

§3. The neutrino emissivity and the specific heat 

Now we are in a position to discuss the cooling behavior by using the above 
results @)'|). First let us calculate the specific heat of the quark matter. To this end, 
we must derive entropy of the quark matter, which is written as 

S = -(l^) =-8E[n p logn p + (l-n p )log(l-n p )]. (15) 

Here we have omitted the contribution of anti-quarks for the sake of brevity. They 
will be taken into account in the numerical calculations in the next section. From 
the entropy, we get the specific heat as usual way: 

*- T (£)„„■ (16) 

Note that another specific heat C p is nearly equal to the Cy in the Fermi gas as is well 
known in statistical physics. The entropy does not contain the derivative of the gap 
variable by virtue of the gap equation (14). Hence it is continuous function of T across 
the critical temperature. On the other hand, the specific heat contains the derivative 
of the gap variable so that it is discontinuous as a function of the temperature, which 
is well known in electron superconductors. The sudden enhancement of the specific 
heat comes from the formation of Cooper pairs. Because a piece of heat applied 
to the superconductor is expended on the destruction of Cooper pairs instead of 
thermal motion. 

Next let us consider the neutrino emi^ivity. It was calculated in the quark 
matter in the normal phase long years ago.cP The result is given by 

C Q m = const.T 6 , (17) 

The temperature dependence of the emissivity can be understood easily. Let us 
investigate the quark URCA process (1). There are u,d quarks and electrons which 
are approximately degenerate with Fermi momentum pp(i). Each Fermion gives one 
power of T from the phase space integral {d?pi — > Pp(i) 2 dEi cx T). Thus we have 
T 3 from the u,d quarks and electrons. In addition, the phase space integral for the 
neutrino gives d 3 p u oc E 2 dE u oc T 3 . One power of T from the emitted neutrino 
energy, E u , cancels a factor T _1 from the energy-conserving 5 function. Since the 
momenta of the degenerate particles are restricted to lie close to their respective 
Fermi surfaces, the angular integrals give no temperature dependence. Altogether, 
we thus have C m oc T 6 . The constant on the right-hand side of the above equation 
depends on the more detailed calculation. Since the value is not important for our 
later discussions, we do not consider it. 
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This derivation of Eq.(17) is based on the assumption that each Fermi energies 
are larger than the temperature, Ep 3> T. Noting that pp oc p 1//3 , this condition 
will be approximately satisfied for the high density quark matter near the critical 
temperature considered in this paper. 

Our interest is now modification of the emissivity due to the color-superconductivity 
and its influence on the cooling behavior. Again see the weak decay (1). There are 
two kinds of d quarks: the (Cooper) paired particles and the unpaired ones. For the 
former case, the d quark might decay in the following mode: 

(ud) — > u + u + e~ + i7 e . (18) 

Noting that Aq ~ lOOMeV < T c , this process is prohibited by the energy conserva- 
tion law. Therefore the /3-decay is suppressed by the formation of the Cooper pairs. 
The rate is estimated to be 

the number of unpaired d— quark 
total number of d— quark 

Similarly the second process of (1), 



(19) 



{ud) + e" -> d + d + u e , (20) 

is suppressed in the same manner as the first one. Consequently the neutrino emis- 
sivity in super phase is given by 



C em = r 2 Cfl (21) 

According to the usual BCS theory &, the rate of the unpaired d quark is given by 
the following equation: 

r = l-4^r^dk. (22) 
3 TkpJo 8E k v ; 

The first term comes from the contribution of the gapless d quarks and the sec- 
ond from those of the unpaired u,d quarks. As a result, the neutrino emissivity is 
suppressed by the pairing correlation. 



§4. Numerical results 



Numerical calculation has been carried out in the case of the chemical potential 
/x = 500MeV. If the coupling constant g is taken as lGeV, the resulting gap energy 
becomes about lOOMeV, which is a typical case used in other references. It is as- 
sumed that the emissivity (21) derived in the previous section is valid under these 
circumstances. 

First the behavior of the specific heat is drawn in Fig.l. The solid line represents 
the result in the super phase and dotted line denotes that in the normal one. We 
can see a discontinuity at the critical temperature representing the phase transition 
of second order. Below the critical temperature, the specific heat decreases with 
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Fig. 1. The specific heat of quark matter as a 
function of temperature. The solid (dot- 
ted) line represents the result in the super 
(normal) phase. 




Fig. 2. The rate of unpaired ci-quarks as a 
function of temperature. The rate is de- 
fined by the equation (22). 



an exponential curve. On the other hand, it is linear in the normal phase so that 
the cooling rate in super phase becomes larger than that in normal phase at lower 
temperature. 

The rate of unpaired particles defined by Eq.(22) is shown in Fig.2. Since the 
Cooper pairs are formed in the super phase, the number of the unpaired quark dimin- 
ishes below the critical temperature. Substituting the specific heat and emissivity 
considered above into Eq.(2), the cooling rate is calculated and the result is drawn 
in Fig.3 as a function of T. We can see the pairing effects on the cooling rate as 
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Fig. 3. The cooling rate —dT/dt — C m (T) /Cv(T) as a function of temperature. 

follows. First the sudden suppression of the cooling rate occurs across the critical 
temperature. This suppression is strengthened by the emissivity. Below the critical 
temperature, the decrease of the emissivity causes overall suppression of the cooling 
rate, in particular at lower temperature. 

Finally some comments are added as follows. The sudden suppression of the cool- 
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ing rate may give a manifest signal for the color-superconductivity in quark matter 
if we could determine the temperature in the core of the stars in any way. However 
it is very difficult task at present. Another possible signal is direct measurement of 
the decrease of the neutrino emissivity, such as measurement of the neutrino energy 
from a supernova. But this would be also very difficult. Anyway the observation of 
the signal is a future problem. 

Second comment is about the other cooling process that gluons excited in hot 
quark matter decay into neutrino-antineutrino pairs: tP 



gluon -^q + q^ZQ^u + u. (23) 

where q + q means a particle-hole state in the quark matter. This process seems to 
give large contribution to the cooling. In the super phase, however, the particle-hole 
state should be replaced by the two-quasi-particle state whose excitation energy is 
larger than 2Aq. Noting that the gluon energy is smaller than T c < 2Aq, this process 
will be prohibited by the energy conservation law. Therefore this cooling process is 
suppressed by the formation of Cooper pairs and may be neglected in the discussions 
on the cooling behavior. 

In conclusion, the neutrino emissivity decreases on account of the color super- 
conductivity and causes suppression of the cooling rate. This phase transition leads 
to a sudden discontinuous suppression of the cooling rate in cooperation with the 
specific heat. 
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